Zero-bias conductance peak and Josephson effect in graphene-NbTiN junctions 
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We report electronic transport measurements of graphene contacted by NbTiN electrodes, which 
at low temperature remain superconducting up to at least 11 Tesla. In devices with a single super- 
conducting contact, we find a more than twofold enhancement of the conductance at zero bias, which 
we interpret in terms of reflectionless tunneling. In devices with two superconducting contacts, we 
observe the Josephson effect, bipolar supercurrents and Fraunhofer patterns. 
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I. INTRODUCTION 

Graphene has attracted attention since its isolation 
from graphite^ and much of its electronic properties 
are now well understood^. Graphene differs from con- 
ventional two-dimensional electron systems by its lin- 
ear dispersion relation at low energies and a vanishing 
band gap. There has been a lot of interest in study- 
ing quantum transport in graphene (G) contacted by 
superconductors (S). So far, most experiments have fo- 
cused on electrical transport in SGS Josephson junc- 
tions. Tunable Josephson supercurrents^, multiple An- 
dreev reflections^^, Andreev bounds states'", phase dif- 
fusion phenomena^, macroscopic quantum tunnelling^, 
and superconducting phase transitions^^' have been 
observed in devices employing Al^'^'^'^^^^, W^, Ta^^, 
p]^7,8,i9^ Pbln^, Sn^^'^^ as superconductors. An inter- 
esting and experimentally yet unexplored direction is to 
demonstrate specular Andreev reflection^^ by realizing 
a superconducting gap larger than the potential fluctu- 
ations in the graphene. Another interesting possibility 
is to study the interplay between superconductivity and 
the quantum Hall effect ^^^^ which requires contacting 
high-mobility graphene with a superconductor with a 
large critical magnetic field. While we were finalizing 
the manuscript, measurements on graphene with Nb and 
ReW contacts at high magnetic field were posted, see 
Ref.26. 

In this study we report electronic transport through 
NbTiN based SGN and SGS junctions, where N is a 
normal, non-superconducting metal electrode. NbTiN 
has a large gap, a high critical temperature (Tc) and a 
high upper critical perpendicular magnetic field {B^^). 
We present electrical measurements in a field-effect ge- 
ometry as a function of source-drain bias, temperature 
(50 mK-5 K) and external magnetic field (0-11 T) ap- 
plied perpendicular to the sample. In SGN devices at 
sub-Kelvin temperatures and up to moderate magnetic 
fields, we observe a zero-bias conductance peak. We an- 
alyze this peak in terms of reflectionless tunneling. ^^^^ 
In SGS devices, we observe gate-tunable supercurrents 
and discuss their magnetic field response, which exhibits 



characteristic Fraunhofer patterns. 

The remainder of the manuscript is organized as fol- 
lows. In Sec. II we describe the details of the sample 
fabrication. We discuss the measurements in Sec. Ill, 
and we conclude in Sec. IV. 



II. SAMPLE FABRICATION 

Graphene samples were prepared on Si substrates with 
a 285 nm thick oxide layer by mechanical exfoliation from 
natural graphite (NGS NaturGraphit GmbH). Graphene 
monolayers were selected by optical contrast.^^ The elec- 
trodes were defined by standard electron beam lithogra- 
phy. The normal contacts consist of a 8 nm Ti adhesion 
layer with a 50 nm AuPd alloy on top. The ^ 30 nm thick 
NbTiN superconducting contacts were made by DC sput- 
tering of a NbTi target (30% Ti, 70% Nb atomic percent- 
age) in a Ar/N2 plasma in a Nordiko-2000 system using 
an unbalanced parallel plate configuration.^^ The Ar and 
N2 flows were 100 seem and 4 seem and the deposition 
pressure was about 6 mTorr. The deposition conditions 
were optimized for producing high quality NbTiN films 
with low stress^^'^^ and a Tc of about 13 K corresponding 
to a BCS superconducting gap of 2 meV. The electrical 
measurements indicate an upper critical perpendicular 
magnetic field Bq^ in excess of 11 T at 50 mK (the exact 
value is not known). 

We investigated three ways of contacting graphene 
with NbTiN: (1) direct sputtering of NbTiN on the 
graphene, (2) sputtering of NbTiN on a Ti protective 
layer and (3) sputtering of NbTiN on a Ti/Au protective 
layer. Direct sputtering of NbTiN on graphene leads to 
very high contact resistances of hundreds of kl] (measure- 
ments of these devices are not discussed further). This is 
attributed to damage to the graphene layer underneath 
the S electrode due to its exposure to the sputtering 
plasma and/or the highly energetic particles involved in 
the sputtering process (DC voltages are of the order of 
380 V). In order to prevent this problem we fabricated 
devices in which we protected the graphene by two differ- 
ent approaches. For the type A (Fig. la) devices, we first 
covered graphene with a thin layer of Ti (10 nm) evap- 
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FIG. 1: Type A SGN device with 7 nm Ti. (a) Schematic 
showing dimensions and electrodes configuration, (b) False 
color SEM image (the white scale bar represents 500 nm). (c) 
dV/dl versus Vgate at Idc = 10 /iA and Iac = 100 nA in two 
gate sweep directions. 



orated in an Eva 450 e-beam evaporator, which involves 
particles with energies of the order of only 1 eV. Next, 
in 3 to 5 minutes the sample was transferred through air 
into the sputtering system, which unavoidably leads to 
oxidation of the Ti. Prior to the sputter deposition of 
NbTiN, the oxidized Ti was cleaned by an Ar RF plasma 
of 200 W for 5 min in a pressure of 6 mTorr. This cleaning 
procedure leaves about 7 nm Ti. For the type B devices 
(Fig. 5a) we e-beam evaporated Ti(2 nm)/Au(2.5 nm) 
as a protective layer. Since no oxide is expected to form 
during the sample transfer, no RF plasma cleaning was 
done before depositing NbTiN. 



III. MEASUREMENTS 

The electrical transport measurements were performed 
in a ^He/^He dilution refrigerator (Leiden Cryogenics) 
with a base temperature of 50 mK using standard low 
frequency lock-in technique which allowed simultaneous 
measurement of the current-voltage characteristic and 
the differential resistance/conductance. The wiring be- 
tween electronics and samples involved pi filters at room 
temperature and RC and Cu-powder filters at the mix- 
ing chamber stage. All measurements were done in a 
two-terminal configuration and the data presented here 
were corrected for the RC filter series resistance of 2.5 kQ 
per filter. 

In Fig. la and lb we show a schematic and a false-color 
scanning electron micrograph (SEM) of a type A device. 
The graphene flake is 2.6 /im wide. The S and N elec- 
trodes completely overlap the flake and are separated by 
about 250 nm (edge to edge). In Fig. 2a we present two- 
terminal differential resistance {dV/dl) measurement at 
50 mK as a function of the back-gate voltage Vgate foi" 
the as-fabricated SGN junction. The measurement was 
done with an AC current of 100 nA superimposed on 
a DC current of 10 /iA. The dV/dl is rather high over 
the entire gate voltage range: given the device dimen- 
sions, the doping level and an estimated carrier mobility 



of ~ 2000 cm^V~^s~^ observed in a nearby device on the 
same flake, the graphene differential resistance should be 
well below 1 kQ. Other devices made in the same batch 
show similar or even higher differential resistances, but 
type A devices with a thicker Ti protective layer (about 
20 nm) and type B devices (without RF plasma clean- 
ing) do not. Three-terminal measurements of the S and 
N contact resistances (to the graphene covered by the 
contacts) show values of 400 Q and 700 ft, respectively. 
Therefore, we conclude that most of the resistance is due 
the transition region from the covered graphene to the 
uncovered graphene (denoted as G'). We believe that 
stress in the NbTiN films and/or RF plasma cleaning 
may cause damage to the graphene area around the S 
contacts. 

In Fig. 2a we show the dl/dV for Vgate = V as a 
function of voltage bias at 60 mK (for two perpendicular 
magnetic field values of T and 11 T) and at 4 K for T. 
We see that dl/dV is non-constant over the entire bias- 
voltage range of ±20 mV. Since the NbTiN gap is only 
^ 2 mV, the non-constant dl/dV above 2 mV cannot 
be attributed to the superconducting contacts. Indeed, 
when normalized to a high bias value of 10 — 20 mV, the 
three dl/dV traces in Fig. 2a fall on top of each other 
except for a bias range of less than 2 mV. The observed 
behavior is also incompatible with known mechanisms 
that lead to a bias-dependent graphene resistance. ^^^^ 
Possibly it is caused by the damaged graphene near the 
contacts. 

Fig. 2b shows dl/dV^s at low bias for several gate volt- 
ages. We observe a zero bias conductance peak for en- 
ergies Ec smaller than 0.1 meV and a conductance dip 
at energies Ed of around 0.2 meV. Both features are 
present at all gate voltages. The zero bias conductance 
is enhanced by more than a factor of two compared with 
its value at ~ 1 mV for almost all the back gate volt- 
ages investigated, see Fig 2c. Ignoring possible magnetic 
moments^^'^^, to the best of our knowledge only reflec- 
tionless tunneling can explain an enhancement a factor 
of more than two. Hereby the transmission of a tunnel 
barrier between a superconductor and a normal metal is 

enhanced due to the diffusive transport in the normal 
metal27 29 

Following the semiclassical approach of van Wees et al. 
we sketch the principle behind reflectionless tunneling 
in Fig. 3. The quasi-particles move from the right nor- 
mal metal reservoir towards the S electrode through the 
diffusive graphene. A potential barrier exists at the 
graphene/ S interface, which in our case is assumed to 
be due to the damaged graphene around the S electrode. 
An electron (ei in Fig. 3) that hits the barrier can be ei- 
ther Andreev reflected as a hole (hi) or normal reflected 
as an electron, continuing along path 2. Due to scattering 
on impurities the electron moving on path 2 has a chance 
to hit the barrier once again where it can retroreflect as 
a hole (h2). Retracing path 2, this second hole reaches 
the initial point where it can undergo normal reflection. 
Constructive interference between the first and the sec- 
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FIG. 3: Illustration of reflectionless tunneling. Quasi- 
particles move from the right reservoir towards the super- 
conductor through diffusive graphene. A potential barrier ex- 
ists between graphene and the superconductor. The Andreev 
reflection probability at an otherwise poorly transparent in- 
terface is enhanced due to interference effects (see text for 
details) . 



FIG. 2: Measurements of the SGN device shown in Fig. 1. 
(a) dl/dV as a function of Vtias for Vgate = V at 4 K 
(S = T) and 50 mK (for B = T and B = 11 T). (b) 
dl/dV vs. bias voltage at various gate voltages indicated 
by labels (T = 50 mK, B = T). The curves are offset in 
40 /iS steps; the top curve has no offset, (c) Zero-bias dl/dV 
(Go) normalized to the value at 1 mV vs. Vgate- The zero 
bias conductance is enhanced by a factor larger than two at 
almost all gate voltages. 



ond hole increases the total Andreev reflection probabil- 
ity of the incoming electron ei. In a diffusive sample 
there are a multitude of such paths and their respective 
contributions add up leading to an enhancement of the 
conductance. At finite bias, the phase of h.2 averages out 
and the conductance enhancement disappears. 

The order of magnitude of the cutoff energy Ec for 
which coherence is lost is determined by the phase co- 
herence time T(f) according to Eq ^ ^/t^. Taking 
Ec ^ 0.065 meV from the dl/dV at Vgate = V, we 
estimate ~ 10 ps which is comparable with the val- 
ues obtained in Ref.^^. The phase breaking length is 
given by = ^JDt^. Here, D is the diffusion constant 
given by D = where / is the mean free path and 

vp ~ 10^ m s~^ the Fermi velocity. Taking = 10 ps 
and a mean free path of ~ 17 nm (corresponding to a 
carrier mobility of 1000 cm^/Vs at 30 V from the Dirac 
point) gives ^ 280 nm, which is comparable to the 
sample length L = 250 nm. 

Now we turn our attention to the magnetic field and 
the temperature dependence. In Fig. 4a we show the 
dl/dV^s as a function of bias voltage for various mag- 
netic fields applied perpendicular to the sample surface. 
A finite magnetic field breaks the time reversal symme- 
try and introduces a phase difference between the inter- 
fering holes. Coherence is lost when the loop formed by 
path 2 and the superconductor encloses one flux quantum 
$0 = h/e. Taking 5 = 4 T as the field where reflection- 
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FIG. 4: (a) dl/dV vs. bias voltage at Vgate = V for 
different magnetic fields, (b) Temperature dependence of the 
dl/dV at Vgate = V and 5 = T. For clarity, the data in 
both panels is offset vertically by incrementally subtracting a 
value of 40 /iS. The topmost curves have no offset. 

less tunneling is suppressed (see Fig. 4), we estimate an 
effective area of 1 x 10^ nm^. At high magnetic fields one 
expects Landau levels to develop as the graphene enters 
the quantum Hall regime. However since we did not see 
clear signs of quantum Hall plateaus, it is likely that the 
disorder in the sample was too high for the quantum Hall 
effect to develop. 

In Fig. 4b we plotted the temperature dependence of 
the dl/dV at Vgate = V and 5 = T. The measure- 
ments were taken after a mild current annealing step^^ 
performed at base temperature. For about 10 minutes, 
we slowly ramped up a DC current applied between the 
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FIG. 5: Type B SGS device with Ti/Au/NbTiN electrodes 
made without RF cleaning prior to depositing NbTiN. (a) 
Schematic of the device structure with the relevant dimen- 
sions, (b) False color SEM image (scale bar 200 nm). (c) 
dV/dl at Vgate = V, B = T. (d) Color plot of dV/dl vs. 
hias and Vgate at B = T. (e) Color plot of dV/dl vs. B and 

^biaS ^9ate = V. 



before deposition of NbTiN. Fig. 5a and Fig. 5b show 
the schematic and a false color SEM image of the device. 
The graphene flake is 1.5 jam wide and the electrode sep- 
aration is 150 nm (edge to edge). In Fig. 5c we show 



the dV/dl as a function of bias current at Vg 



gate 



V 



and B = T. A critical current Ic of 4 nA is observed. 
We note that Ic is relatively small given the junction 
dimensions and the carrier density. Also, we did not ob- 
serve supercurrents in type B junctions with 280 nm or 
larger electrode separations, possibly due to the weak- 
ening of the proximity effect in the Ti/Au bilayer. The 
sharp peaks in dV/ dl we interpret as self-induced Shapiro 
steps since the energies involved closely match those of 
standing waves formed in the metal box enclosing the 
sample, which has a size of 6.0 cm.'^^'^'^ In Fig. 5d we 
show dV/dPs as a function of the gate voltage and bias 
current as 2D color plot. The dark region corresponds to 
a supercurrent through graphene. The magnitude of the 
critical current depends on the charge carrier density and 
decreases when going from metallic conduction towards 
the charge neutrality point (at ~ 22 V), in agreement 
with previous reports^. In Fig. 5e we show dV/dls as a 
function a perpendicular magnetic field and bias current 
at zero gate voltage. The critical current is modulated by 
the magnetic field revealing the well-known Fraunhofer 
pattern. The area extracted from the periodicity of the 
Fraunhofer pattern amounts to 0.16 /im^ in close agree- 
ment with the geometrical area of 0.2 /im^ deduced from 
the SEM image. The width of the first lobe is smaller 
than twice the period of the higher order lobes which 
indicates non-uniform current flow. 



S and N contacts, up to a current density of 4.5 A cm~^. 
While this lead to a ~ 50% increase in the overall conduc- 
tance, the zero-bias conductance peak and the broader 
conductance dip remained at about the same energies. 
This behavior is similar to what was reported in Ref.^^. 
From Fig. 4b we observe that the zero-bias peak drops at 
about 1 K. This thermal energy scale is consistent with 
Ec extracted from the width of the peak. 

Now, we briefly discuss measurements of SGS junc- 
tions. In a type A SGS device with similar dimensions 
as the SGN device discussed earlier and made in the 
same batch, no supercurrent flowing through graphene 
(Josephson effect) was observed. We understand this in 
terms of the poor transparency of the G-G' interfaces 
which hinder the diffusion of the Cooper pairs into un- 
covered graphene. In another batch, with a thicker Ti 
layer (about 20 nm after the RF plasma cleaning), we 
did observe bipolar supercurrents and Fraunhofer pat- 
terns in several graphene junctions for electrode spacings 
of up to 400 nm^^. This is indicative of less damage to 
the G-G' transition region. These data are not discussed 
further. 

In Fig. 5 we show measurements of a type B device. 
Given the Ti(2 nm)/Au(2.5 nm) protective layer, no 
oxide is expected to form during the sample transfer 
through air so we skipped the RF plasma cleaning step 



IV. CONCLUSION 

In this paper we reported electrical transport measure- 
ments of SGN and SGS junctions with NbTiN supercon- 
ducting electrodes. In SGN devices with a barrier region 
near the S contacts, we observe a zero-bias conductance 
peak, which we think may arise from reflectionless tun- 
neling. While we are unable to identify its origin un- 
ambiguously, this peak has a striking resemblance to the 
recently measured zero-bias anomaly in a candidate topo- 
logical superconductor Cua^Bi2Se3.^^ In junctions with 
transparent interfaces, supercurrents were observed for 
electrode separations of up to 400 nm. By improving the 
fabrication procedure we believe that high carrier mo- 
bility graphene samples and clean graphene-NbTiN in- 
terfaces can be obtained for experimentally investigating 
specular Andreev reflection and the interplay between 
Andreev reflection and the quantum Hall effect. 
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